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1. Introduction 



The present work is yet another chapter in the continuing story, inspired by the 
seminal papers of the interrelation between quantum integrability and combinatorics. 
It concerns the combinatorial interpretation of the ground state components of a quantum 
integrable ID Hamiltonian, or equivalently of the equilibrium distribution of a Markov 
process on planar connectivities of points (also called link patterns). We refer the reader 
to [Q for an overview of the relevant model, related to the Temperley-Lieb alegbra. These 
ground state entries turn out to be integers which count classes of Fully Packed Loop 
(FPL) configurations, better known in the mathematical literature as Alternating Sign 
Matrices (ASMs). An important detail is that while the first articles considered 
models with periodic boundary conditions, subsequent work also studied also types of 
boundaries conditions ]|,|5]||[7| . In particular in |3|,[7] one can find a definition of the so- 
called mixed boundary conditions which will be used here. In the loop model to ASM 
dictionary, boundary conditions of the model translate into symmetry classes of ASMs. 

Here we make use of the method of study developed in [|J, and of its extension con- 
sidered in 0, or more precisely of its reformulation in terms of qKZ equation given in [1C]. 



The essence of the method is to generalize the problem by introducing extra parameters 
(in a similar spirit as ||11|| ), making the model inhomogeneous. We go further by replacing 
the eigenvector equation with a more general linear system, which we call a bit loosely 
the qKZ equation, which lets us vary an extra parameter q. The latter generalization, al- 
though not strictly necessary for our purposes, is both conceptually satisfactory (allowing 
to distinguish clearly what properties are dependent or independent of q), and of some 
mathematical interest, in relation to representation theory |J and algebraic geometry |T(J . 

The plan of the article is as follows: in Sect. 2 we give a brief description of the loop 
model with mixed boundary conditions and define our notations. In Sect. 3 we define 
the qKZ equation and discuss the relevant solution. Sect. 4 contains the main application 
of this formalism: the computation of the sum rule of the model with mixed boundary 



conditions, as well as its refinement in the sense of [12"|. Sect. 5 provides concluding 
comments. 
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2. Definition of the model 



We provide here the minimum about the model with mixed boundary conditions 
introduced in that is needed for our purposes, as well as other useful definitions. Let 
L be a fixed positive integer. Consider points on a line, numbered from 1 to L, plus 
an extra point to the right of them. The space of states consists of right- extended link 
patterns (or right-extended matchings, in the language of H), that is pairings of the 
points {1, . . . , L} between themselves, allowing for an arbitrary number of unpaired points 
which are connected to the extra rightmost point, in such a way that all these connections 
can be drawn in the half-plane without crossings, see Fig. 1. There are (^ 2 j) sucn nn k 
patterns. 



L = 3: 




1 2 3 




1 2 3 



1 2 3 



L = 4: 



12 3 4 




12 3 4 




12 3 4 



12 3 4 



12 3 4 



12 3 4 



Fig. 1: The right-extended link patterns of size L = 3,4. 



These link patterns form the canonical basis of a vector space, on which the one- 



boundary-extended Temperley-Lieb algebra [|f3||14|1 acts. Its generators, ei,...e^_i and 
/ send link patterns to other link patterns according to the standard graphical rule ex- 
emplified on Fig. 2. Each time a loop is formed, it must be erased and the state must 
be multiplied by a weight r. On the other hand "boundary loops" that go through the 
rightmost point produce no weight. The operators and / satisfy the following relations 
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Fig. 2: Examples of action of Temperley-Lieb operators. 



The term "mixed boundary conditions" refers to the fact that link patterns are closed 
on the left, but open/extended on the right. 

Set r = —q — q~ x . One can define the i?-matrix 



Ri(z) 



(qz-l/q)I+(z-l)ei 
q-z/q 



(2.2) 



where / is the identity operator, and the boundary i?-matrix [15,16 



K{z) 



( Z - q yC)(z-C/q)I+(l-q)(l-z 2 )f 
(qz-(/q)(z-q/C) 



(2.3) 



where £ is an additional free parameter of the model (this expression can also be found in 
a different parameterization in the appendix of 0). 
These operators satisfy the unitarity equation: 



Ri(z)Ri(l/z) = I K(z)K(l/z)=I 
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(2.4) 



the Yang-Baxter equation: 

Ri(z)Ri +1 (zw)Ri(w) = R i+1 (w)Ri(zw)R i+1 (z) (2.5) 

and the boundary Yang-Baxter (or reflection) equation: 

R L _ 1 (w/z)K(z)R L _ 1 (l/(wz))K(w) = K(w)R L _ 1 (l/(wz))K(z)R L _ 1 (w/z) (2.6) 

as well as the usual commutation relations [Ri(z), Rj(z')] = for \i — j\ > 1 and 
[K(z),Ri(z')} = 0forz<L-l. 

An integrable model is usually defined by a infinite family of commuting transfer 
matrices. In the present context it is more relevant to consider scattering matrices Si, 
% = 1, . . . L — 1, defined by 

Si(z u ... , z L ) =Ri{s Zi/z i+1 )R i+1 (s Zi/z i+2 ) . . . Rl-i{s Zi/z L ) 
K{l/{ SZi )) 

R L -i(sz L Zi) . . . Ri(s z i+1 Zi)Ri-i(s Zi-xZi) . . .Ri(sz\Zi) 
R 1 (z l /z 1 )...R t - 1 (z t /z l - 1 ) (2.7) 

where s is yet another free parameter. Intuitively (as in the Coordinate Bethe Ansatz), Si 
corresponds to a "particle" starting at location i with spectral parameter Zi, moving to the 
left and crossing other particles with spectral parameters Zj, j < i, then reflecting at the 
left boundary by having its spectral parameter changed to l/(szi) (but with a trivial left 
boundary i?-matrix), then going back to the right crossing all other particles, reflecting to 
the right with the matrix K and having its spectral parameter inverted to s Zi, and finally 
moving back left to its original location by crossing particles with spectral parameters zj, 
j > i. 

It is a simple exercise to show, using the properties above, that these operators Si 
satisfy the following commutation relations: 

Si(zi, . . . , s zj, . . . , z L )Sj(z 1 , . . . , z L ) = Sj(zi, . . .,szi,.. . , z L )Si(zi, . . . , z L ) (2.8) 

In particular note that when s = 1, the Si commute. This is the "physical" situation 
that leads to the diagonalization problem for the Si (which will be investigated in detail 
in Sec. 4). However, for the time being, we do want to keep s free. 
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3. gKZ equation and its Laurent polynomial solution 

3.1. The linear system 

Consider the following set of equations for a Laurent polynomial in the variables 
J with values in the vecter spaee spanned by the right-extended link patterns of 

size L: 

Ri(z i+1 /zi)^ L (z u ...,Zi, z i+1 , ...,z L ) = ^l(zi, • • • , z i+ i, Zi,..., z L ) (3.1a) 
K(z L )^ L ( Zl , ...,z L ) = V L (z u z L -u (3.16) 

Vl{zi,...,Z L ) = ^ L (l/(sZ 1 ),Z 2 ,...,Z L ) (3.1c) 

These equations mimic the definition of the model, corresponding to the behavior of 
the ground state eigenvector f l at the bulk, the right boundary and the left boundary 
respectively. In order to explain more precisely the connection with the previous section, 
one can compute S^l for a solution of equations (|3.1] ). Going back to the definition 
(P-7|) of the operators Si, we apply successively Eqs. ( |3.1a| ) for z,i + 1, . . . , L — 1, then 



Eq. ( CT) , then Eqs. (|3TTaQ for L- 1, . . . , 1, then Eq. (gTj) , then Eqs. ( g7Ta| ) for 1, . . . , i-1 
and the end result is 

Si(zi,...,Z L )^ L (zi,...,Z L ) = ^ L (zi,...,szi,...,z L ) (3.2) 

Thus, $i is a vector that has very simple transformation properties under the action of 
the Si. Eq. ( |3.2| ) is in fact the gKZ equation (or more precisely, a variant of it since the 
original gKZ equation [lTTJ is really the analogue for periodic boundary equations). 



3.2. Laurent polynomial solution 
We now claim the following: 

If s = q 3 , there exists a solution of Eqs. ( |3.1| ), which is a centered Laurent polynomial 
of total degree width 2L(L — 1)/2 and degree width in each variable 2(L — 1) . The solution 
is unique up to normalization and can be constructed explicitly by starting from the base 
component (empty link pattern with no pairings, first in the examples of Fig. 1 ) 

*L-,o = J] - q~ X ){q~ l Zi - qs-h- 1 ) (3.3) 

l<i<j<L 



1 \I/ l can be made into a polynomial by multiplying it by (zx ■ • ■ ZL) k for k sufficiently large; 
this however makes Eqs. ( |3.16, c ) look slightly less nice. 
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and applying Eqs. ( |3.1a, Uj) . 



Similar statements have been made in various related models, and the strategy to 
prove them always follows the same general pattern. Here we shall simply explain why the 
equality s = q 3 , as well as Eq. ( |3.3| ), necessarily hold, and how the latter fixes all other 
components. 

First we rewrite Eqs. ( |3.1a, fl ) by using the explicit forms (Eqs. (|2.2|), (^73])) of R and 
K. We obtain the following dichotomies: 

al) If the link pattern n contains the pairing (z, i + 1), then noting that eiix = ttx we can 
rewrite the component 7r of Eq. ( |3.1a|) as 



*tt' = (QZi ~Q 1 z l+1 )d i ^ n (3.4) 

7r'^7r:ei7r'=7r 

where d\ is the divided difference operator: diF(zi, Zi+i) = Ei^l+li^ll — _^j Zl1 z%+1 ^ . 
a2) If i and i + 1 are not paired in 7r, then \E f 7r U i+1=g 2 2 . = and ^ n /(qzi — q~ 1 Zi + i) is 
symmetric in Zi, Zi+±. 
and: 

bl) If the point L is unpaired in tt, then noting that fn = nwe can rewrite the component 
7T of Eq. (|3TT5D as 

£ *fV = ^(9 - g -1 ^ 1 )^ - ?C _1 ) (3-5) 

where dF(^) = F ^ - F ^ . 
b2) If the point L is paired in n (i.e. not connected to the rightmost point), then fy n has 
zeroes at Zl = q~ 2 (, zl = (?C _1 - an d ^-n-/((q — Q~ 1 C z l 1 )( z l ~qC~ 1 )) 1S unchanged 
by z L -> 

Using recursively these properties, as well as the remaining equation ( |3.1c| ), one can 
determine all the "trivial" factors in the components depending on the properties of 7r: 

★ If no pairings occur between the consecutive points {*,... , j}, Oi<fc<z<j (*? z k~q~ l z i) I ^n- 
■k If no pairings occur between the consecutive points {1, . . . , j}, YliKkKKji^ 2 ~ s z k z i) I ^n- 

★ If L is unpaired and there are no pairings between the consecutive points {«,..., L — 1} 
and no unpaired points among them, rii<fc<z<L(l ~ Q 2z k z i) I ^tt- 

★ If L is paired and there are no pairings in the interval {*,..., L}, Ylk=i( ( i~ ( i~ 1 C z k 1 )( z k ~ 
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The first and second properties imply that ^l-o contains all the factors of Eq. ( |3.3|) , 
and since these exhaust the degree width of the claim, there can be no more (monomials 
included, since ty no must be centered). Furthermore, it is clear that Eqs. ( p.4j) and ( |3~5l ) 
allow to build all other components starting from ^l-o, and can only preserve or lower 
the degree width. Finally, we shall show explicitly in the following examples how s = q 3 
is necessary for L = 2, L = 3. Below we shall exhibit recurrence relations that allow to 
produce ^l-2 from so that the condition s = q 3 must be true for any L. 

EXAMPLE: In size L = 2, there are only two components: the empty link pattern (num- 
bered 0) and the link pattern that pairs 1 and 2 (numbered 1). Starting from 

*2;0 = {qziz^ 1 - q~ 1 )(q~ 1 Z2 - qs^z^ 1 ) 
and noting that / sends 1 to 0, we compute 

* - 1 (n n-^t7- 1 \(7 at- x \m - (9 4 - s)(gz 2 - c/g)(z 2 - g/Q 

V 2;l - [q-q <^Z 2 ){Z2-q(, JCW 2; o - 7- r-g 

1-q (l-q)q 2 sz 2 
But we also have that e\ sends to 1 so that we can also compute 

n _vD ( n7 n~ x 7\c)^ - ( g -Q 3 )(Q 2 ^ -Z2){\-q 2 z x z 2 ) 

U — ^2;0 — \qZ\ — q Z 2 JClW 2) l — 7- r— 5 

(1- q)q 6 sz 1 z 2 

so that we conclude that s = g 3 . Substituting this leads to 

^2;0 = {Z IZ2 ~ 1 - q~ 2 ){z 2 - g" 1 ^- 1 ) 

*2;l = -q~ 2 (q - <rW)Cz2 -qC' 1 ) 

(note in particular the simplification in ^ 2 -± which becomes a Laurent polynomial in q). 

In size L = 3, we use the ordering of Fig. 1. Starting from \I>3 ; o, we compute ^ 3 -i = 
jz^iq-q^Cz^izs-qC'^d^s-fi and # 3;2 = -#3,0 + (<? z 2 -q~ 1 z 3 )d 2 '^ 3 . 1 ; skipping these 
intermediate steps, we note that since e 1 " 1 (2) = {0, 1, 2}, we must have 

= ^ 3 ;0 + ^3;i-(g^i-g" 1 22)9i*3,2 = (s - g 3 )(g 2 Zl - z 2 )(g 2 Zl - z 3 )(g 2 z 2 - z 3 ) R(z u z 2 , z 3 ) 

where R is some irrelevant rational fraction. We conclude again that s = q 3 , and after 
simplification: 

^3;0 ={ z i z 2 l ~ Q~ 2 )( z 2 ~ Q~ 1 Zi 1 )(z 1 z^ 1 - q~ 2 )(z 3 - X ) (z 2 z 3 ~ 1 - q~ 2 )(z 3 - q~ 1 z 2 ~ 1 ) 

*3;i ={zi z 2 l ~ Q~ 2 )(Z2 ~ q~ l z^){q - q- l Qzz l ){z 3 - qQ~ X )x 

x (-q~ 2 (z 1 + z 2 ) + q- 3 (l + q~ l )z 3 - q'^z' 1 + z' 1 ) + q~ 3 (l + q'^z' 1 ) 

*3;2 = (z 2 Z 3 ~ 1 - q~ 2 )(z 3 - q~ 2 Z^ X ) (-q^zt (z 2 +z 3 )+q- 2 (l+q- 1 )z 2 z 3 +(q- 3 C+^ 1 )((.l+q)zi (22+23)) 

-(i-g~ 1 ) 2 g" 1 (i+g" 1 )-9- 4 2r 1 (2 2 - 1 +23- 1 )+q" 2 (i+q- 1 )2 2 z 3 - 1 ) 
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3.3. Recurrence relations 

We derive here a recurrence relation with respect to the size L satisfied by ^l- 
Consider ^ L\z i+1 =q 2 Zi for 1 < i < L — 1. Using the results of previous section we know 
that only the components with a pairing (i, i + 1) are non-zero. Defining the operation <pi 
that inserts into a link pattern n of size L — 2 the pairing (z, i + 1) (thus shifting the points 
i, . . . , L — 2 by 2), we claim that 

*£,,&(7r)(2lj • • • , ^)|z i+1 =g 2 Zi = PiO^ Zl, • • ■ , • • • , Zi) *Z,-2,7r0l, • • • , • • • , Z L ) 

(3.6) 

where Pj is some Laurent polynomial of degree width 2(2L — 3) which is symmetric in 
Zi, . . . , Zi-i, Zi+2, %l- The common factor Pj can be found by factor exhaustion, up to 
a multiplicative constant; it is then simple to check that the l.h.s. of Eq. ( |3.6j ) satisfies 
Eqs. (|3.1|) at size L — 2; one concludes by uniqueness of the solution of the prescribed 



degree that Eq. ( |0|) holds. Remains a multiplicative constant that needs to be computed 
by taking an example. 

Since we shall need in what follows only one recurrence relation, we do the computation 
of Pi only for i = L — 1. In this case note that the empty link pattern of size L — 2 is sent 
by 4>L-i to the link pattern with one pairing (L — 1, L). The latter is the only antecedent 
of the empty link pattern by /; thus its component is 

*L,l = y^(q - q- l Qzl l ){z L -qC^Lfi 

Setting zl = q 2 ZL-\ we note that ^l,q vanishes, so that we are left with only one term in 
d^ L ,o- Explicitly, 

^L,i\ ZL = q ^ L ^ ={qz L -x ~ r X )(l " q-\zl\) II - 4~ 2 )^ - TV 1 )* 

l<i<j<L-l 

x II ^ z l-i - i)( z l-i - q~ A h x ) 

\<i<L-2 

Dividing by \1/ l-2,o(zi, . . . , zl-2) we obtain 

PL-iizL-^Z!, . . ., z L -2) = (qz L -i - C _1 )(l - ^C^Z-i)* (3.7) 

L-2 

x II ( ZiZ i~-i - q~ 2 )(zL-i - q~ 1 z~ 1 )(z l zl^_ 1 - q)(z L -i - q~ 4 z~ X ) 
i=i 

All the factors can be expected on general grounds (zeroes of \Pz,), but our computation 
has produced the normalization constant. 
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4. Sum rule at the stochastic point 

4-1- Stochastic point 

A value of the parameters that is of special interest is q = e ±2l7r / 3 , that is the weight 
of a loop is r = 1, and s = q 3 = 1. At this special value we notice that all operators e$ and 
/ have the co-vector v = (1, . . . , 1) as left eigenvector, with eigenvalue 1: vei = vf = v. 
So do the operators Ri, K, as well as S{. In stochastic processes, this is interpreted as the 
conservation of probability. Also, we note that the operators Si now commute according 
to Eq. (|2.8|) . It will therefore come as no surprise that Eq. (|3.2| ) becomes the equation for 
their common (right) eigenvector with eigenvalue 1: 

Si(z!, . . . , z L )^ L (zi, ...,zl) = ^l(zi, ...,zl) (4.1) 

We can now consider the limit Zi — > 1, in which all operators Si tend to the identity. 
Expanding to first order in any variable we note that -^-Si\ Zk= \ = aijH + j3ij, where 
and flij are irrelevant constants, and H is the Hamiltonian 

L-l 

H=J2a + af ( 4 - 2 ) 

i=l 

in which we have set a = 3/(£+ 1 + C -1 )- H, being a sum of L operators with the same left 
eigenvector, also possesses it, with eigenvalue L — 1 + a. Furthermore, if a is real positive, 
it is easy to show that the matrix of H in the canonical basis satisfies all properties of 
the Perron-Frobenius theorem, so that L — 1 + a is its largest eigenvalue, and its right 
eigenvector defined by 

H^ L = {L-l + a)^ L (4.3) 

is its Perron-Frobenius eigenvector. Of course $i is up to a multiplicative constant 
. . . , 1), but we choose its normalization in such a way that the component indexed 
by the empty link pattern is a L Z/ / 2 J . 

Let us first discuss the case Q = 1, i.e. a = 1. All componenents turn out to be 
integers, and in is given a remarkable conjectural combinatorial interpretation of these 
numbers (a la Razumov-Stroganov conjecture ||) which we explain now. 

Define a Fully Packed Loop configuration (FPL) to be a coloring of edges of a square 
lattice in two colors (say, occupied and empty) in such a way that external edges are 
alternatingly occupied and empty and that vertices have exactly two occupied adjacent 
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edges and two empty adjacent edges. FPLs are known to be in bijection with Alternating 
Sign Matrices (ASMs), see for details. A FPL is called vertically (resp. horizontally) 
symmetric if it is symmetric with respect to the vertical (resp. horizontal) axis. Now we 
can state the conjecture: with our normalization each component is an integer and 
counts the number of Horizontally and Vertically Symmetric Fully Packed Loop configu- 
rations of size 2L + 3 with the connectivity of the corresponding link pattern it, see Fig. 3. 
As a corollary, the sum of components v ■ \& = Y2 n ^ * s ec l ua l to the total number of 
such configurations, or equivalently the number of Horizontally and Vertically Symmetric 
Alternating Sign Matrices (HVSASMs). 




3 4 

Fig. 3: Example of a HVSFPL and the corresponding link pattern. 

Let us now go back to a general value of a. A refinement of the previous conjecture, 
formulated in [|12], states that the components provide an a-enumeration of the HVSFPLs, 
the power aL fc / 2 J being related to the number k of paths crossing the horizontal symmetry 
axis say, to the left of the center, that go straight down two steps below this axis. Again, 
the sum of components becomes the a-enumeration of HVSASMs, where a weight a is 
given to each pair of — l's on the row below the horizontal symmetry axis and to the left 
of the center. 

Remark: at a = the Hamiltonian becomes that of the model with closed boundary 
conditions (whose ground state eigenvector is related to the counting of VSASMs or even- 
sized analogues, see [|,|6|,|3|] ) , and one expects to recover the results of []T8| , which discusses 
the inhomogeneous generalization. Keeping only the leading term of \& as £ — > 0, one finds 
that only components with all points paired (except one if L is odd) are non-zero, which is 
the sector preserved by the operators ei, i = 1, . . . , L — 1; and one could naively expect these 
limiting components to be exactly the polynomials introduced in |fT8| , up to a common 
factor. Interestingly, it is not the case - if only because the polynomials in our case have 
higher degree and are coprime. A related crucial mismatch is that the gKZ equations are 
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different: that of the closed boundary conditions has a shift s = o 6 , as opposed to s = q 3 
here. Only at q = e ±2l7r / 3 , where the two \& are ground state eigenvectors of the same 
transfer matrix or scattering matrices, is the identification correct. Indeed one can check 
that the components of ^l(C ~ > 0) develop a common symmetric factor at q = e ±2l7F / 3 
(which is precisely the function Xl(zi, ■ ■ ■ , z£) to be defined below, cf Eq. |(4.7)| ). A similar 
phenomenon of interrelation between s = q 3 and s = g 6 solutions of qKZ was observed in 
P,|T9|1 in the case of periodic boundary conditions. 



EXAMPLE: For L = 4, with the order of link patterns as on Fig. 1, 



H 
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1 
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1 

Vo 
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1 





1 

1 + a 

1 






a 
1 

2 





1 


1 + a 
1 



°\ 




1 

a 

2/ 



(4.4) 



The sum of each column is 3 + a, as should be, and ^4 is 

#4 = (a 2 , 3a, 2a(3 + a), 3, 3a(2 + a), 3(2 + a)) 



(4.5) 



The corresponding HVSASMs are drawn in |I2"|] . At a = 1, ^4 = (1,3,8,3,9,9), which 
sums to 33, the total number of HVSASMs of size 11. At a = 0, #4 = 3(0, 0, 0, 1, 0, 2), 
which is the ground state eigenvector of the Hamiltonian with closed boundary condi- 
tions, normalized in such a way that the GCD is the total number of VSASMs of size 5. 
Each component, once divided by this GCD, counts VSASMs of size 5 with prescribed 
connectivity. 

4-2. Sum rule 

We are particularly interested here in the "sum rule" Zl(z±, . . . , zl) '■= v ■ 
Zl) = 5^tt ^L,n(zi, . . . , zl) where we have set q = e ±2l7F / 3 but retained the 
spectral parameter dependence. Zl is a Laurent polynomial of degree width 2L(L — l)/2, 
2(L — 1) in each variable. Applying v on the left to Eqs. (|3.1| ) (remembering that v is 
left eigenvector of Ri, K with eigenvalue 1), we find immediately that Zl is a symmetric 
function of its arguments, and that it is invariant by z% —>■ 1/zi. Next, we sum over n in 
Eq. (|3.6|) and obtain 



Z L (zi, ...,z L ] 



z L =q*z L 



.! =Zl-2(zi, ■ ■ ■ , Zl-2)X 

x (l-qzL-iC 1 )^ 



(4.6) 



L-2 



MIT 1 



qz L -iz i 
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where we have made use of Eq. (3"77). The recurrence relation is started with Zq = Z\ = 1. 



Zl-, being a Laurent polynomial of degree width 2(L — 1) in say zl, is entirely fixed by its 
value at 2L — 1 points. Eq. (|4.6| ) (using the symmetries of Zl) provides us with the values 
of Zl at zl = q^z^ 1 , which is more than enough. Therefore we only need to find the 
unique solution of this recurrence. 

Introduce the following determinant ratio: 

detfz J+rj/21 " 1 -z~ j ~ im+ \<- ■< 

Xn(z!,..., Zn) = ^ - JJ ^ (4-7) 

det(^ - Zi J )i<i,j<n 

This function, a character of the symplectic group Sp(2n), already appears in [TJ] and 
in 1 20|j , as we shall discuss in more detail. It is a symmetric Laurent polynomial of its 
arguments, invariant by Z{ — > 1/-Z£, of degree width 2|~n/2(n/2 — 1)], and it satifies the 
recurrence relation: 

n-2 

Xn(zi, • • • , z n )\ Zn=q 2 Zn _ 1 = Xn-2(zi, z n - 2 ) JJ(1 - q Z n -iZ~ l )(Zi - q 2 z~^ x ) 

1=1 

We then easily find that 

Z L (zi, ...,z L ) = Xl(zi, . . . , z l )xl+i(z!, ...,z l ,() (4.8) 

since it has the right degree and satisfies the recurrence. 
Example: The first few sums: 

Z 2 = z x + z 2 + C + Zi 1 + z^ 1 + C 1 
Z 3 = (zi +z 2 + z 3 + Z1 1 +Z2 1 + z 3 ~ 1 )x 

x [z\z 2 + Z1Z3 + z 2 z 3 + z± 1 z 2 ~ 1 + z± 1 z 3 ~ 1 + z 2 ~ 1 z 3 ~ 1 + (C + C~ 1 )(zi + z 2 + z 3 + z^ 1 + z^ 1 + z^ 1 ) 
+ z 1 z 2 ~ 1 + z 1 z 3 ~ 1 + z 2 z 3 ~ 1 + z 2 z± 1 + Z3Z1 1 + Z3Z2 1 + 3) 



Finally, we can take the homogeneous limit Z{ = 1. In all that follows we rede- 
fine the normalization of Z so that the smallest component is ^0 = a\ L ^ 2 ^ as before: 
Z L (a) = 3" L(L - 1) / 2 aL L / 2 Jz L (l, . . . , 1). First we set ( = 1 (a = 1) as well. The limit where 
parameters go to 1 in Eq. ( |4.7| ) can be found, as always with these types of determinants, 
by setting Z{ = exp(ei) and letting e go to zero. This is strictly equivalent to using Weyl's 
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dimension formula for representations of Sp(2n). Define for future use hi = i + \i/2] — 1; 
then x„(l, ••.,!) = ni<i<j<„ i J=r IL<i<j<„ ttt- so that 

X„(l,...,l)3-r»«»/ 2 -t)l = JJ L^ + lJg)*' = i, 1,2,3,11,26, 170,... 

2fn-l-fc 

where we have divided by the appropriate power of 3 so that the resulting number has the 
following interpretation: it is, for n even, the number of Cyclically Symmetric Transpose 
Complement Plane Partitions on a n x n x n hexagon, and for n odd, the number of 
Vertically Symmetric Alternating Sign Matrices of size n. These numbers already appear 
as the sum of components in the case of closed boundary conditions @JT^]. 

Going back to Z^, the powers of 3 are cancelled by those coming from \&£ ; o(l, . . . , 1) = 
gL(L-i)/2^ anc j we conclude that 

^ A L3fc/2 + lj(3fc)!fc! 

\ = 11 (2fc + l)!(2fc)! = 1j 2 ' 6 ' 33 ' 286 ' 442 °' • • • 
which is the formula for the number of HVSASMs of size 2L + 3 conjectured by Mills 



and proved by Okada [|0 



If we now reintroduce the £ dependence, we find an expression which, as already 
stated, is conjecturally the weighted sum of VHSASMs where a weight a is given to each 
pair of — l's on the row below the horizontal symmetry axis to the left of the center. In 
particular, in |T^] was formulated a conjecture about pl = logZi,| a= i (which should 
be the average density of such —l's). Let us explain here how to compute this quantity. 

One starts with the following identity related to GL(N) characters, valid for arbitrary 
integers hi. 

X h? f n < N - 1 



Hn{ku ■ ■ ■ ' M = g Twt-^) = i «„-* + i(*i, • • • , h N ) n > N — 1 ^ 

where s n is the GL(N) Schur function associated to the Young diagram with one single row 
of length n. This is proved by considering the generating function ^2 H n t n , using the fact 
that H (hi, . . . , hN, = and concluding with the identity J2 n >o s n(hi, . . . , h^t 71 = 
ex P Sn>i(Si h™)t n /n. Here we only need the case n = N, which, according to Weyl's 
dimension formula, corresponds to computing the first logarithmic derivative of a GL(N) 
character wrt one parameter at the identity of the group, the result being which 
is nothing but the number of boxes of the Young diagram. 
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The corresponding "symplectic" identity can essentially be obtained by replacing h{ 
with tij. In the present case we wish to compute the second derivative with respect to £ 
of logXL+i (since the first vanishes). We deduce, with n = N = L + 1, 

d 2 . ...... i ^ 



d<2 log Xi+1 (l, . . . , 1, Ok-, = (L + 1)(2L + 3 ) E ^ - ( 4 ' 10 ' 



where we have specialized to /ij = i + |~i/2] — 1 . 

Performing the summation: J2f=i( h i ~ * 2 ) = ( L + 1 ) L L / 2 J (5[(L + 1)/2J + 2)/3, and 
changing variables to a = 3/(£ + 1 + C _1 ): n °t forgetting the normalization \l/o = a^ 2 -', 
we obtain 



which is equivalent to the expression given in [12 



5. Comments and conclusion 

5.1. Relation to UU ASM partition junctions 

In view of the conjecture of U and its refinement in [12|, it is tempting to try to 
identify our sum rule with parameters Zl(z\, ■ ■ ■ , zl) with the partition function of 
Alternating Sign Matrices with two U-turn boundaries as introduced in [|IT| - that is the 
partition function of the six-vertex model with its usual weights depending on spectral 
parameters (and q = e ±2l7r / 3 ) and certain special boundary conditions. The UUASMs 
are a slightly more general class of ASMs than VHSASMs, but they are the right objects 
to consider once spectral parameters are introduced; one can recover VHSASMs from 
UUASMs in the homogeneous limit, as will be discussed below. The comparison reveals 
subtle differences. UUASMs being defined in even size only, we must discuss the two 
parities separately: 

★ For even L, depends on L spectral parameters, just like ours (in the notations of [|TT 



and [^0( , these are the square roots of our parameters Zi), plus two boundary parameters, 



called b and c. In fl20|l , the case b = c = q 1 ! 2 is treated and the partition function is found to 



be equal (up to irrelevant prefactors) to Xl(zi, • • • , zl)xl+i(zi, ■ ■ ■ > Zl, 1). This coincides 
with Zl at £ = 1. It is in fact not too hard to check that if one sets b 2 = c 2 = q/£, one also 
obtains the dependence on £. The amusing observation is that in the homogeneous limit 
Zi = 1, the condition b = c = q 1/2 selects VHSASMs (of size 2L + 3) among UUASMs; 
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however setting b 2 = c 2 = q/Q does not. This provides an apparently non-trivial equality 
of the weighted sum of UUASMs, where the weights are on the two U-turn boundaries 
(more precisely, each U-turn gets a weight of q — (,/q if it is the correct orientation of the 
arrow for a VHSASM of size 2L + 3, 1 — £ if it is the wrong orientation) and the weighted 
sum of VHSASMs where the weight a = 3/(£ + 1 + C _1 ) is on the row below the horizontal 
symmetry axis (for each pair of — l's to the left of the center). As a check one can consider 
the case ( = q 2 , which selects among UUASMs the VHSASMs one size less (2L + 1): the 
conjectured equality implies that the number of VHSASMs of size 2L + 3 with maximal 
power of a is equal to the number of VHSASMs of size 2L + 1, which can indeed be proven 
(in fact there is a bijective proof, which is left as an exercise). 

* For odd L, the situation is even more mysterious. ZY+i depends on L + 1 parameters, 
that is one more than in our model. In ||20|| , the case b = c = q~ x l 2 (which, in the 



homogeneous limit, selects VHSASMs of size 2{L + 1) + 1 = 2L + 3) is treated and the 
result is that ZY+^zi, . . . , z L +i) oc Xl+i(z\, . . . , zl+i)xl+i(zi, • ■ ■ , z L +i), where xl+i is a 
certain character of 0(2(L + 1)) which satisfies xl+i(^i, ■ ■ • , 1) oc Xl(zi, • ■ ■ , zl). Thus 
Zl+i coincides with Z^{C, = 1) on condition that one of its spectral parameter is set to 
1. It is unclear how to reintroduce ( in the partition function of UUASMs, and inversely 
unclear how to incorporate the extra spectral parameter in our model. 

5.2. Limit q — > 1 and geometry 

The main focus of the present paper is on the special value q = e ±2l7F / 3 . However 
the rational limit of the qKZ equation is also interesting. This corresponds to the limit 
q — > 1 (in [10 1 the limit is indifferently q — > ±1, but this is because the shift s = q 6 ] here 



s = q 3 — » 1 imposes q — > 1), that is expanding to first non-trivial order in h as q = e h t 2 , 



Zi = e hWi . In jTIJ, it was found that solutions of the usual type A (periodic boundary 



conditions) rational gKZ equation based on Hecke algebra quotients are related to the 
geometry (more precisely, equivariant cohomology) of certain sl(L) orbital varieties. This 
idea was generalized in |2^] to Temperley-Lieb models with other boundary conditions 
which are apparently much closer to the present work; the models were shown to be 
related to type B, C, D orbital varieties. However a key difference is that the geometric 
interpretation was ensured from the start by the fact that the .R-matrix (both bulk and 
boundary) dependence in the spectral parameter had a specific form which matched the 
Hotta construction of the Joseph representation for orbital varieties. This resulted in 
slightly unusual variants of the Temperley-Lieb algebra. Here the strategy is different: 
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using the standard 1-boundary extended Temperley-Lieb algebra, we are forced to use a 
certain boundary .R-matrix whose geometric interpretation is not obvious. Note that the 
example of [p3l , |24j1 shows that more exotic varieties/geometric constructions than orbital 
varieties/Hotta construction can appear in integrable models. We hope to come back to 
this question soon. 

5.3. Positivity conjectures 

In [II]] it was observed that in the case of periodic boundary conditions, the compo- 



nents of solutions of the qKZ equation in which the variables Z{ are set to 1 are polynomials 
with positive coefficients in the variable: (1) r = —q — q~ l in the case s = q 6 ; and (2) 
t' = q 1 ! 2 + q~ 1 / 2 in the case s = q 3 . See also p5[ for an interesting conjecture on the 



combinatorial interpretation of the sum of homogeneous components for arbitrary q in the 
case s = q 6 . 

We list here the first few homogeneous = C = 1) in terms of r' and after dividing 

by some common factors: 

#2 = (M) 

* 3 = (l,2,2 + r' 2 ) 

# 4 = (1, 3, 2(3 + r /2 ), 2 + r' 2 , 5 + 3r /2 + r' 4 , 7 + 2r /2 ) 

Clearly the positivity still seems to hold. Note that the special values considered above, 
namely q = e ±2l7F / 3 and q = 1, correspond to r' = 1 and r' = 2 respectively. 

5.4- Other boundary conditions and prospects 

Several times in this article we have mentioned the similarities with the work IT? 



This is no surprise because the present model is a generalization of the one considered 
in 



P~8|| : the Hamiltonian of the model with mixed boundary conditions has an additional 
boundary term, and when the coefficient a in front of it is set to zero one recovers the 
Hamiltonian of the model with closed boundary conditions. At the technical level, the 
main difference is that the boundary term implies that one must include a non-trivial 
boundary i?-matrix satisfying the reflection equation. (Another difference, irrelevant at 
the "physical" value q = e ±2l7F / 3 , is that the shift s in the qKZ equation turned out to be 
q 3 here, whereas it is q G in ||18|| ). 

Similarly, one could add another boundary term on the left, which results in the model 
with open boundary conditions 0, an even broader generalization. The analysis of the 



corresponding qKZ equation is more involved, and is deferred to a future publication |26[ . 
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